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In the paper we study the nonlinear mechanical cooling processes in the intrinsic quadratically 
optomechanical coupling system without linearizing the optomechanical interaction. We apply the 
scattering theory to calculate the transition rates between different mechanical Fock states with the 
use of the resolvent of Hamiltonian, which allows for a direct identification of the underlying physical 
processes, where only even-phonon transitions are permitted and odd-phonon transitions are forbid¬ 
den. We verify the feasibility of the approach by comparing the steady-state mean phonon number 
obtained from transition rates with the simulation of the full quantum master equation, and also 
discuss the analytical results in the weak coupling limit that coincide with two-phonon mechanical 
cooling processes. Furthermore, to evaluate the statistical properties of steady mechanical state, we 
respectively apply the Mandel Q parameter to show that the oscillator can be in a nonclassical me¬ 
chanical states, and the phonon number fluctuations F to display that the even-phonon transitions 
favor to suppress the phonon number fluctuations compared to the linear coupling optomechanical 
system. 

PACS numbers: 42.50.Wk, 42.65.-k, 07.10.Cm 


I. INSTRUCTION 

Realization of cooling of mechanical motion within the 
vicinity of quantum regime is a key ingredient for the 
broad range of its applications, including quantum infor¬ 
mation processing (QIP) [1-3], high-precision measure¬ 
ment [4, 5], and probe of quantum behaviour of macro¬ 
scopic system [6], etc. Recently, theoretical and experi¬ 
mental studies are mainly focused on the linearized treat¬ 
ment of linear coupling optomechanical system, where 
the coupling term is proportional to the mechanical dis¬ 
placement X [7-9]. However, there are some litera¬ 
tures have begun to generalize the investigations on the 
nonlinear mechanical cooling processes, such as second- 
sideband laser cooling of trapped ions [10] and cool¬ 
ing in the single-photon strong optomechanical coupling 
regime [11], which can lead to nonthermal steady states, 
and even the nonclassical, sub-Poissonian states, and 
quantum systems far from thermal equilibrium hold great 
promise for the investigation of fundamental physics and 
the implementation of practical devices [12, 13]. 

By now a new type of optomechanical system, i.e., 
quadratic optomechanics, has been proposed, where the 
coupling term is proportional to square of mechanical 
displacement x^, which is realizable in the setups of 
membrane-in-the-middle optomechanical system, ultra¬ 
cold atomic ensembles and superconducting electrical cir¬ 
cuit systems [14-16]. Distinct from the displacement 
of mechanical equilibrium position in the linear cou¬ 
pling optomechanics, the quadratic coupling will modify 
the resonant frequency of the mechanical oscillator [17- 
19]. Thus the quadratic optomechanical system can be 
used for quantum non-demolition (QND) measurement 
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of individual quantum jumps [20]. In addition, the ro¬ 
bust stationary mechanical squeezing and electromag- 
netically induced transparency (EIT) from two-phonon 
processes in quadratic optomechanical system are also 
achievable [21, 22]. For the practical applications one 
should minimize the influence of thermal noise, and lin¬ 
earized approach on the cooling of quadratically coupled 
mechanical oscillator has been studied in Ref. [23] . With 
the advancement of experimental technology, now it will 
become attractive to theoretically investigate the nonlin¬ 
ear cooling processes induced by the intrinsic quadrat¬ 
ical coupling. For example, utilizing avoided cross¬ 
ings effects in a multi-mode optical cavity containing 
a flexible dielectric membrane enables the significantly 
enhanced quadratical coupling strength that can reach 
> 30MHz/nm^ [24, 25]; in the superconducting elec¬ 
trical circuit quadratic optomechanical system, the cou¬ 
pling strength g compared to dissipation rate k could be 
raised to g/K > 0.1 [16]; cavity optomechanical system 
with ultracold atomic ensembles also provides an obvious 
characterization of quadratic-coupling regime [15]. These 
platforms offer the possibility to enter the single-photon 
strong quadratically optomechanical coupling regime. 

Here, we investigate the nonlinear mechanical cooling 
processes in the intrinsic quadratically optomechanical 
coupling regime without linearizing the optomechanical 
interaction. We apply the scattering theory to calculate 
the transition rates between different mechanical Fock 
states with use of the method of resolvent of Hamilto¬ 
nian [26, 27]. This approach allows for a clear identifi¬ 
cation of the underlying physical processes, where only 
even-phonon transitions are permitted while the odd- 
phonon transitions are prohibited. In the paper we first 
derive the mechanical transition rates, verify the consis¬ 
tence of transition rates with the results of the simulation 
of full quantum master equation under the consideration 
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of thermal mechanical clamping, and then discuss the an¬ 
alytical results in the weak coupling limit that coincide 
with two-phonon mechanical cooling processes [28, 29]. 
To indicate the difference to the linear coupling optome¬ 
chanical system by evaluating the Mandel Q parameter, 
it tends to achieve the nonclassical mechanical states in 
quadratical optomechanics. Furthermore, we apply the 
phonon number fluctuations F to present that the even- 
phonon transitions favor to suppress the phonon number 
fluctuations. 

The paper is organized as follows. In Sec. II the 
quadratical optomechanical system is introduced, in 
Sec. Ill the scattering theory on the single-photon 
quadratically optoemchanical cooling process is pre¬ 
sented and the cooling limits in the weak and strong 
coupling regimes are discussed, and in Sec. IV the non- 
thermal and nonclassical statistical properties of the me¬ 
chanical oscillator is analyzed. At last the conclusions 
are given. 


II. DESCRIPTION OF PHYSICAL MODEL 

We consider the quadratically coupled optomechanical 
system in which an optical cavity mode is parametrically 
coupled to the square of the position of a mechanical 
oscillator, and the Hamiltonian is written as 


pressions of coupling constant can be referred in lit¬ 
eratures [30, 31]. The interaction gives the losses of the 
cavity at rate k, with the definition 


K = 



2 

S{ujk - ujr). 


( 3 ) 


The mechanical dissipation rate is usually much smaller 
than the cavity damping rate, and we will take it into 
consideration in the following. 


III. RESOLVENT METHOD ON NONLINEAR 
COOLING PROCESSES 

In this paper we will resort to the scattering theory to 
explicitly evaluate the transition rates which quantita¬ 
tively determine the cooling and heating processes. This 
approach will allow for a direct identification of the un¬ 
derlying physical processes. For convenience we separate 
the total system Hamiltonian according to H = Hq + V, 
where the redehned main operator Hq combines Eq. (1) 
and free electromagnetic field, i.e., 

Hq = -Aa^a + uJmb'^b + ga^a{P + b^ + i?emf, (4) 
and the small interaction part is 


Hq = a + uJrny^b + ga^ aiV^+ bY, (I) 


V = W, + H', 


( 5 ) 


where lor is the frequency of the cavity mode d, LOm is 
the frequency of the mechanical mode 6, and g is the 
quadratically optomechanical coupling strength between 
the cavity mode and mechanical motion of the oscillator. 
In order to make the present system stable, the Routh- 
Hurwitz criterion implies that the condition uj^+Asg > 0 
should be satisfied, in which s is the number of photons 
inside the cavity. In the frame rotating at the driving 
laser’s frequency lur, the drive is described by the Hamil¬ 
tonian H' = VL{a^ + a) with the driving strength fl. We 
are looking forward to the regime H ^ k, where k is 
the damping rate of the cavity field, and thus the cavity 
states with more than one photon can be neglected. 

The damping of the cavity field is described by the 
interaction between the mode a and the environmental 
mode Cke, where Hem! = - WL)cj],gCfce is the free 

dynamics of the electromagnetic field external to res¬ 
onator in the reference frame rotating at the laser fre¬ 
quency ujR. Here, the sum runs over all modes, identified 
by the wave vector k and orthogonal polarization e, the 
operators Cke and annihilate and create photons in 
the corresponding mode, and ujk = ck denotes the mode 
frequency. The cavity-environment interaction 

^ g^k} -b acl^) (2) 

ke 

accounts for the coupling of the cavity mode with the 
modes of the external radiation field and the explicit ex¬ 


in which A = lur — ujR is the detuning between laser ujr 
and resonator frequency lor. In the case of weak optical 
drive strength H, it is sufficient to work to second order 
of H in the transition rate to well describe the scattering 
processes. Here we suppose that the initial state is \i) = 
jO; n; 0^^), where the cavity mode is in the vacuum state 
jo), jn) is phononic excitations and jOfce) is environmental 
vacuum mode. We are interested in the processes that 
change the motional state of the membrane induced by 
one photon followed by dissipation to environment, and 
denote the final states ]/) = lO;TO;lfee). Thus we only 
need to retain the transition amplitude proportional to 
H, which is given by 

Tfi = {f\VGQ{E,)V\i), (6) 


where Go{z) = 1/{z—Hq^) is the resolvent of the effective 
Hamiltonian Hq^ = Hq — under the consideration 

of cavity damping, and Ei = nuim is the energy of initial 
state [32]. With the help of polaron transformation 


U = exp 




( 7 ) 


and the relation 

- Hf )U^] = I, ( 8 ) 

Z Hq 
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where I is the identity matrix, the resolvent of the effec¬ 
tive Hamiltonian can be reexpressed as 


Goiz) = 


z- Hf 


= 


z-Hf 




= U^ 


U{z - Hf)W 


u. 


( 9 ) 


By applying the form of polaron transformation, the 
transformation of the effective Hamiltonian becomes 


U{z 


= z — — lS.a)a + ujm\l^ +'^a^ag/wr, 


X 1/2) 


1 

2^,71 - *-a'a 


( 10 ) 


With the explicit form of interaction part V in Eq. (5) 
and using the initial state |i) and final state |/), we can 
find out that the nonzero part of transition amplitude is 

r» = (/I |i) = 

= {n+^)ujm -{1+ ^)w£\ 

= Wmx/TTdg/W/, (11) 


in which Si^n = (^|exp —^(JP — b^^)/2 |n) with f = 

j h\{l+4,g/LOm) is the coefficient of squeezed number state 
in Fock representation and defined as [17, 18, 33] 


Floor[|] Floor[f ] _ 

(coshO”+i/2 2^ k'\k\ 

—0 A /—0 

X (12) 

where the function Floor [a;] gives the greatest integer less 
than or equal to x, and is the one-photon coupled 
membrane’s resonant frequency [19]. 

Then the corresponding transition rate is defined as 


i?fi = 27r^|rfi|"(5(F7-F/), (13) 

fee 


where i5-function guarantees the energy conservation and 
the sum covers all relevant polarizations and wave vec¬ 
tors. The transition rate becomes 


(k) 

9k, 


Rn =27rfl^ 

ke 

X S(ujk - ujr) 


E 


^1,171^1,- 

Wy(i^ny+A+2I 


(14) 


With the definition of cavity damping in Eq. (3), the 
transition rate which denotes the change of the motional 


state of mechanical oscillator from n to m phonons be¬ 
comes 


77/2 — ^ n^m^n — 






^ ( 5 (i)(/,n) + A-f if 


(15) 


Equation (15) presents a clear physical view on the pro¬ 
cess of an incident photon scattered by the quadratic 
cavity-optomechanical system with changing the motion 
of mechanical oscillator from n to m phonons. When a 
laser photon turns into the cavity photon, the phonon 
state changes from the initial number n to an inter¬ 
mediate number I demanding that \l — n\ is even due 
to the term of two-phonon operators, and the ampli¬ 
tude is inversely proportional to the frequency detuning 
, n) -|- A which relates to the photon-coupled mem¬ 
brane’s resonant frequency ujm\ and proportional to the 
matrix element of squeezed number state As the 

photon leaves the cavity and dissipates into the contin¬ 
uum of modes in free space, it induces a transition in 
mechanical motion from I to m phonons which also de¬ 
mands that \l — m\ is even, and its amplitude is indepen¬ 
dent of driving detuning but dependent on the cavity- 
reservoir coupling strength and also determined by 
the coefficient of squeezed number state Si^m for the two- 
phonon term. The output photon possesses the energy 
tOL + [n — m)uJm and carries away (n — m) phonons, in 
which (n — m) should be even since \l — n\ and |Z — to| 
are both even. 

Now we stress the differences of nonlinear cooling rates 
between quadratical and linear coupling optomechanical 
systems. The quadratical optomechanical coupling mod¬ 
ifies the mechanical frequency of the oscillator while the 
linear coupling just displaces the mechanical equilibrium 
position. Thus the frequency detuning of the transition 
is closely dependent on the coupling strength 
g while the effective detuning in the linear coupling op- 
tomechanics is independent of the coupling strength in 
the single-photon level. Moreover, the transition am¬ 
plitude of the linear coupling optomechanics is propor¬ 
tional to the Franck-Condon overlap factor for the single¬ 
phonon operators while the transition amplitude is pro¬ 
portional to the matrix elements of squeezed number 
state for the two-phonon operators. Therefore, the form 
of transition rate in Eq. (15) is different from that of lin¬ 
ear coupling optomechanical system in Ref. [11, 27], and 
may be helpful to generate more interesting mechanical 
states. 

Taking into account of the thermal damping of me¬ 
chanical oscillator with rate 7 m and thermal phonon 
number nth, the set of rate equations for the mechani¬ 
cal oscillator becomes 


Pn=- 7 mnth(n + P)Pn “ 7 m(nth + l)nPn 

+ 7 mnthnP„_i -I- 7m(nth + l)(n -|- l)P„+i 
— ^ r„_tmPn + ^ hm-triPm, (I m — n| is even) 

m^n m^n 

(16) 
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where Pn is the phonon number distribution in the n 
phonon state. 


A. Cooling process in the resolved-sideband and 
weak quadratical coupling limits 


The occurrence of cooling processes in the strong 
quadratical coupling regime can absorb multiple two- 
phonons, and modify the eigenfrequency of mechanical 
oscillator simultaneously, which can lead to the genera¬ 
tion of uncommon mechanical states. To understand this 
behavior we first consider the case of resolved-sideband 
limit ujrn ^ K and weak quadratical coupling limit 
g <C ojm , where the scattering is mainly dominated by the 
two-phonon transitions r„_>„± 2 , and the rates in Eq. (15) 
are explicitly simplified as r „_,.„_2 = n{n — l)r^ and 
r „_,.„_ i_2 = (n -b l)(n -|- 2 )r.^ in these limits, in which 


r = 

(k;/2)2 + (A + 2u;™) 2’ 
^ («:/2)2 + (A - 2u;™)2 


(17) 


describe the strengths of two-phonon absorption and 
emission processes. The results are consistent with ef¬ 
fective two-phonon quadratically mechanical cooling the¬ 
ory [28], in which the cavity photon is linearized and 
then served as the two-phonon reservoir in the weak op¬ 
tomechanical coupling limit. The set of rate equation in 
Eq. (16) becomes 


Pn =7m(nth + 1)[(?^ + l)^’n-|-l “ nP„] 

7m^th[(^ “b l)Tn TlP^—l] 

+ r^[(n -b l)(n -b 2 )P „+2 - n(n - 1)P„] 

- rt[(n + l)(n + 2)P„ - n(n - 1 )P„_ 2 ]. (18) 


For the high-Q mechanical oscillator and low initial 
temperature, i.e., ^ ignore the ther¬ 

mal damping effects. The cooling process turns to be 
the pure two-phonon transitions, and the phonon number 
distribution can be analytically solved and expressed 
in the form 


P2l\j = (1 - r)r-il+J - J = 0,1, (19) 


where r = T^/T^, and 7 characterizes the relative weight 
of the odd phonon states determined by the initial con¬ 
ditions [34]. In the resolved-sideband regime, when the 
detuning satisfies A = —2uJm, the two-phonon absorp¬ 
tion process dominates the emission process, i.e., r <C 1 , 
and the phonon number distribution is concentrated at 
zero- and one-phonon states Pg = 1 — 7 , Pi = 7 , 
which indicates that two-phonon cooling processes pre¬ 
serve the phonon-number parity, leading to the initial 
odd phonon states cooled to the one-phonon state and 
the even phonon states cooled to the zero-phonon state 
independently. 


In the strong two-phonon absorption regime T^ 3> 
'ynth 3 > r^, which is achievable with A = —2uJm in the 
resolved-sideband limit, finally we can obtain [34] 


Po = 
Pi = 


1 + 2^ 

l-b3C 


l-b3C 


+ 0{^), 
+ 0 (^), 




nth 

«th + 1' 


( 20 ) 


Thus the minimal mean phonon number (n) = 1/(4 -b 
1/nth), and the Mandel Q factor that is defined as Q = 

- 1 with (h) = Ylu=g'nPr,,{n^) = J2u=o'^'^Pn 
becomes Q « —{n), which shows sub-Possionian phonon 
distribution. Thus the membrane is in a nonclassical me¬ 
chanical state. 


B. Cooling process in the strong coupling limit 

In the strong coupling regime higher-order mechanical 
sidebands appear and the multiple two-phonon transi¬ 
tions begin to work. In Fig. 1 we plot the steady-state 
phonon number Uss = (bPb) as a function of detuning 
A = Wi — LOji with different coupling strengths gjujm- 
We observe that in the weak coupling and the resolved- 
sideband limits, i.e., gjuJrn = O.I, the optimal cooling 
occurs at the detuning A = —2uJrn and the cooling limit 
reaches the value close to 0.25, which is in accordance 
to Eq. (20). In the strong coupling limit, the efficient 
cooling processes occur at the multiple two-phonon res¬ 
onances A = {l,n), where the multiple dips of the 

evolution of rigs indicate the occurrence of the multiple 
two-phonon transitions, which will further influence the 
statistical properties of the mechanical oscillator. 

To indicate the validity of the set of rate equations 
obtained by scattering theory in Eq. (16), we compare the 
steady-state phonon number with numerically solutions 
of the quantum master equation 

+ ^^[a]/3-b7m(nth + l)C[b]p 
+'rmnthCp]p, (21) 

performed with the Quantum Optics Toolbox [35] , where 

the Hamiltonian is H = —Aa"^a + ijJm,b^b + ga^a{b + P)^ + 
0(d -b a'^), C[o\p = opo^ — \{o^d,p} is the Lindblad op¬ 
erator of photon (phonon) dissipation. In Fig. 1, we plot 
the results of master equation (21) in variant style of 
lines, which are respectively red solid, blue dashed and 
green dash-dotted lines, for different coupling strengths 
glwm = 0.8,0.4,0.1, and specific markers, which are 
green squares, yellow circles and magenta diamonds, are 
obtained by solving the set of rate equations (16). These 
two results are well matched, the rate equation approach 
coincides with the master equation, and thus scattering 
theory is adoptable in describing nonlinear quadratical 
cooling processes. 
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FIG. 1. (Color online) Steady-state phonon number 
Uss = (b^b) in single-photon quadratical optomechanical sys¬ 
tem as a function of detuning A with the coupling strengths 
5/w m — 0.8 (red solid line and green squares), — 0.4 

(blue dashed line and yellow circles) and gjiiJm = 0.1 (green 
dash-dotted line and magenta diamonds). The other pa¬ 
rameters are = 10~®, nth = 10, = 0.25 and 

Q./k = 0.4. (In the figure the lines are obtained from the 
simulation of the master equation (21) by using the Quantum 
Optics Toolbox and the markers are obtained by solving the 
rate equations in Eq. (16).) 


IV. STATISTICAL PROPERTIES OF STEADY 
MECHANICAL STATE 

The occurrence of the multiple two-phonon transitions 
in cooling processes, especially in the strong coupling 
regime, decreases the membrane’s vibrating energy con¬ 
taining several two-phonons simultaneously, and together 
with the modification on resonant frequency of membrane 
due to the quadratical coupling, the mechanical oscilla¬ 
tor can present different statistical properties compared 
with linear coupling and linearized optomechanics. In 
the section we mainly resort to Mandel Q parameter and 
phonon number fluctuations F to discuss the mechanical 
statistical properties. 

By using Mandel Q parameter it is convenient to char¬ 
acterize nonclassical states with negative values, which 
indicate a sub-Poissonian statistics and have non clas¬ 
sical analog. In Fig. 2, we plot the Mandel Q param¬ 
eter as a function of detuning with different quadrati¬ 
cal coupling strengths. In the weak coupling regime, i.e. 
g/iOm = 0.1, the minimum value of Q becomes negative 
that is close to the value — (h), which agrees with the re¬ 
sults of Eq. (20). With the increasing coupling strength, 
Q can present the negative value as well. For example, 
at gjum = 0.4 it takes a smaller value. Thus the me¬ 
chanical oscillator is in a nonclassical state. However, in 
the strong linear coupling optomechanics [11, 27], it is 
relatively difficult to present the nonclassical properties 
in weak driving regime. 

In the linearized treatment of mechanical cooling 
regime, the rate equation satisfies the detailed balance 
condition which implies that the final phonon state is 
in a thermal state. For the thermal state, (h) = nth, 



FIG. 2. (Golor online) Mandel Q parameter in single-photon 
quadratical optomechanical system as a function of detuning 
A with the coupling strengths g/LJrn = 0.8 (red solid line), 
gjuim = 0.4 (blue dashed line) and g/LJrn = 0.1 (green dash- 
dotted line). The other parameters are 'ymlijJm = 10“®, nth = 
10, KjijJm = 0.25 and H/k = 0.4. 

(h^) = 2njj, -h nth, yielding the phonon number fluc¬ 
tuations F = {b^b^bb)/{b^by = 2 [36]. Here the rate 
equations (16) obviously do not obey the detailed bal¬ 
ance condition and the steady mechanical state will be 
in a nonthermal state. It is more clearly exhibited in the 
single two-phonon cooling process of the weak coupling 
limit in Eq. (18). In Fig. 3 we plot the phonon number 
fluctuations F as functions of detuning A with different 
coupling strengths. Since phonon number fluctuations 
F = 2 for the thermal mechanical state, in the nonlin¬ 
ear intrinsic quadratical optomechanics we can both de¬ 
crease and increase the number fluctuations {F < 2 and 
F > 2), and the fluctuations can be suppressed for a large 
range of detuning. Moreover, compared with the optome¬ 
chanical cooling induced by single-photon linear coupling 
where the fluctuation is mainly enlarged while the me¬ 
chanical cooling occurs here with the mainly decreased 
phonon number fluctuations, which possesses potential 
applications. The suppression is induced by the (mul¬ 
tiple) two-phonon cooling processes that decrease the 
phonon occupation in the zero- and one-phonon number 
states, and compared to one-phonon cooling processes it 
is more achievable to suppress the phonon number fluc¬ 
tuations [11, 28]. 

V. CONCLUSIONS 

In conclusion, we have studied the nonlinear mechan¬ 
ical cooling processes in the intrinsic quadratically op¬ 
tomechanical coupling regime without linearizing the op¬ 
tomechanical interaction. We apply the scattering theory 
to calculate the transition rates between different me¬ 
chanical Fock states with the use of the method of re¬ 
solvent of Hamiltonian, since the approach can present 
a direct identification of the underlying physical pro¬ 
cesses. Due to the preservation of phonon number par¬ 
ity, only even-phonon transitions are permitted and the 
odd-phonon transitions are forbidden. In the paper we 
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FIG. 3. (Color online) Phonon number fluctuations F in 
single-photon quadratical optomechanical system as functions 
of detuning A with the coupling strengths gjixim = 0.8 (red 
solid line), g/ijJm = 0.4 (blue dashed line) and gjijJm ~ 0.1 
(green dash-dotted line). The other parameters are = 

10“®, nth = 10, = 0.25 and fl/n = 0.4. 


first derive the phononic transition rates, and verify the 
feasibility of scattering approach by comparing with the 
results of simulation of the full quantum master equation. 
We also discuss the analytical mechanical cooling limits 
in the weak coupling limit, and find that they coincide 
with two-phonon mechanical cooling processes. Finally, 
the statistical properties of the mechanical state is pre¬ 
sented. The mechanical state can be in a sub-Poissonian 
distribution that characterizes the nonclassical state, and 
the even-phonon transitions favor to suppress the phonon 
number fluctuations for a large range of detuning as well. 
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